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EFFECT OF BODY PERTURBATIONS ON HYPERSONIC FLOW OVER SLENDER 

POWER LAW BODIES 

By Harold Mirels and Philip R. Thornton 


SUMMARY 

Hypersonic-slender-body theory , in the limit as the 
free -stream Mach number becomes infinite , is used to 
find the effect of slightly perturbing the surface of 
slender two-dimensional and axi symmetric power 
law bodies. The body perturbations are assumed to 
have a power law variation (with streamwise distance 
downstream of the nose of the body). 

The perturbation equations formulated herein can 
be used for a variety of problems . Tn particular , the 
effect of boundary-layer development, very small 
angles of attack , and nose blunting can be found . 
Numerical results are presented for (1) the effect of 
boundary-layer development on the flow over two- 
dimensional and axi symmetric slender power law 
bodies, (2) the effect of very small angles of attack 
(on two-dimensional power law bodies ) ; and (3) the 
effect of blunting the nose of very slender wedges and 
cones . 

Differential equations for finding the effect of a 
power law lateral perturbation of the centerline of a 
slender power law body are formulated . No numeri- 
cal results are given . Probably, the most important 
application of these equations is to determine the 
How about axially symmetric power law bodies at 
very small angles of attack. 

INTRODUCTION 

Tn viscid hypersonic flow over slender power law 
bodies was studied in references 1 (o 3. These 
references assume 6 2 C 1, 1/(A/<5) 2 C1 where M 
is the frec-stream Mach number and 5 is a char- 
acteristic streamline slope. Tn the limit 1 f(M&) 2 = 
0, the shock shape and body shape are similar and 
t lie equations of motion can be reduced to a set of 
ordinary differential equations. Numerical solu- 
tions of those “zero-order” equations are tabulated 


in references 1 and 2. Approximate analytical 
solutions arc derived in reference 1. The first- 
order effect- of small but nonvanishing values of 
\/(Md) 2 is also found in references 1 and 2 by ex- 
panding the equations of motion in terms of 

Tn the present paper attention is restricted to 
the limiting ease l/(A/5) 2 — 0. Two-dimensional 
and axisymmetric body shapes of the form 1 
r 6 ~,r m +er v+w are considered where r b is the body 
ordinate, x is distance from the nose (in the frec- 
stream direction), m, N , and e arc constants, and e 
is small. For e— 0, the body is of the simple 
power law type considered in references 1 to 3. 
The flow corresponding to 0 is termed the 
“zero-order” flow herein and may be found from 
the zero-order solutions presented in references I 
and 2. For e small, but not zero, the additional 
term represents a small power law perturba- 

tion of the zero-order body shape. The solution 
for the resulting flow field perturbations is termed 
the “perturbation solution” and is the subject of 
the present- report. The equations developed 
herein ran he used to find the effect of houndary- 
layer development, very small angles of attack, 
and nose blunting on slender power law bodies at 
hypersonic speeds. Numerical results are pre- 
sented for a variety of eases. 

Tt should he recalled (e.g., ref. 1) that the 
assumptions incorporated in hypersonic slender 
body theory are violated in the immediate vicinity 
of the nose (jr— 0) of slender power law bodies 
(except for wedges and (‘ones). Hence, the solu- 
tions found herein are not expected to be valid at 
the nose. However, they are expect ed to he valid 
downstream of this region. 

1 A more specific representation of tlu* Body shapes considered herein is 
given by ( qs. (13) and (G3a). 
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ANALYSIS 

The equations of motion for hypersonic flow 
over slender bodies are utilized to obtain the 
zero-order and perturbation equations for hyper- 
sonic flow over slightly perturbed power law 
bodies. Expressions for shock shape, pressure dis- 
tribution, and drag are then noted. Finally, an 
analytic solution of the perturbation equations is 
given for one class of perturbations. A numerical 
solution of the perturbation equations is required, 
in general. Symbols are listed in appendix A. 
Many of the details of the analysis are relegated 
to appendixes B to IT. 

HYPERSONIC SLENDER BODY THEORY 

The equations of motion for hypersonic flow 
over slender bodies (e.g., ref. 4) are summarized 
in this section. The equations are applicable pro- 
vided that 5 2 <d and l/(A/6) <0(1). The sum- 
mary is the same as that, given in reference 1 and 
is repea ted here for convenience. 

Dimensional variables are barred herein {u, v } 
x } 7, etc.) . See figure 1 for some of these quantities. 



Ficutve 1. -Physical quantities for study of hypersonic 
flow over slightly perturbed power law bodies. Subscript 
zero refers to unperturbed power law body. Note, 
5— /?o(L) / Ij. 

Let 5 represent a characteristic body or streamline 
slope and L represent a characteristic strcamwisc 
length. Two-dimensional and axisymmetrie flows 
are considered, with (7,7) and (u,v) being the 
stream wise and transverse coordinates and veloc- 
ities, respectively. In order to obtain the hyper- 
sonic slender bod}” equations of motion, the follow- 


ing nondimonsional quantities are introduced 
(following ref. 4) : 

x=x/L u~ (u — p = piyM 2 5 2 J) <M 

r = rJL5 v=~v/u m & 

The body shape and shock shape are denoted by 
r b =r b (x) and 7? — 7?(J), respectively, so that 

r>=rj7j U=Ti/L5 ( 2 ) 

If these quantities are introduced into the equa- 
tions of motion, and terms of order <5 2 are neglected 
(compared with 1), the hypersonic slender body 
equations are obtained. These are (ref. 4): 


Continuity: 


dp ,dpr P l n 

-J— _ per — == U 

dr dr r 

(3a) 

/-Moment mil : 


'(c+'£> t-o 

(3b) 

Energy: 


aC/'V) , l .Mp!p ')_ 0 
dr dr 

(3c) 

The boundary conditions are: 
At body surface: 


dr 6 

(4a) 

I T pstream of shock: 


u m = v „ = 0 

(4b) 

Pa> — 1 / (yA/ 2 5 2 ) 

(4c) 

P»= 1 

(4d) 

Downstream side* of shock: 



(4c) 



] (4g) 


Here — 0 for two-dimensional flows and <r= 1 for 
axisymmetrie* flows. This system of equations 
can be solved independently of the ^-momentum 
equation, and therefore the latter is neglected. 
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In the present report it is further assumed that 
1 / (A/’S ) 2 — >0 so that equations ( 4 c) and ( 4 e) to ( 4 g) 
become 


The present problem may be viewed as that of 
finding the flow field associated with shock shapes 
of the form 


0 

II 

8 

( 5 a) 

/^x"(l + to 2 a JV *) 

( 9 ) 

2 d R 

( 5 b) 

where e, a> } and A T are constants and e 

is small. 

ls 7+1 (1j 

To effect a solution, new independent 

variables 



are introduced according to tl 

u* relations 


2 /d 7 ?V 

( 5 c) 




Pa =u±r\ 

( 5 d) 

r r 


(10) 

Ps \r— 1/ 





FLOW ABOUT SLIGHTLY PERTURBED POWER LAW BODIES SO that 


The hypersonic slender body equations of l-lu 1 - 
previous section are now used to find the hyper- 
sonic flow over slightly perturbed slender power 
law bodies. The limiting case 1 / (A/5) 2 -->0 is as- 
sumed. Xote that cither the assumption that 
< 5 Z <0 or that l/(d/ 5) 2 ->0 is violated at the nose 
(# — 0) of a slender power law body (except for 
wedges or cones) so that the solutions found herein 
are not expected to be valid in the immediate 
vicinity of the nose. However, they are expected 
to be valid downstream of this region. 

The zero-order body shape and shock shape ran 
be expressed as (from ref. 1) 

r b ,o(x) = Ti bi oCl' m (Ca) 

7? 0 (5)=CI* (6b) 


where i7 M , C, and m are constants. For the 
remainder of the report the characteristic length 
L is taken to be the stream wise length of the body, 
while the characteristic slope is taken to he the 
zero-order shock ordinate at x— L divided by /., or 


£ _ Rq(L) 7*6,0 (Z>) 

L Lrj bf o 


( 7 ) 


Equations (C) can then be written, in nondimen- 
sional form, 

r bt o=vi>M m (8a) 

II 0 =x m (8b) 

The zero -order flow field is the same as that of 
references 1 and 2 and is considered known. The 
constant ?? M , which is the ratio of the zero-order 
body ordinate to the zero-order shock ordinate, is 
denoted by the symbol in reference 1 and is 
tabulated therein for various values of 7, <7, and m. 


d d mrj d 

d7""5f T^v 

jd = J_ _d 
dr dr; 


( 11 ) 


Those are the same independent variables used 
in reference 1 . The zero-order shock location 
corresponds to rj = 1 while the zero-order body 
location corresponds to = In the new no- 
tation, equation (9) becomes 


+ (12a) 

The boundary conditions (oqs. (5b) to (5d)) 
suggest the following forms for the dependent 
variables: 


i' — m!~ m Vi) 

(12b) 

p m 2 ^ 3(m ~ 1) (Fo+epTj 

(12c) 

P = <Pu\-t$ X 'p2 

(12d) 


where <p, F, and \j/ are functions of 7] and the sub- 
script 0 indicates the zero-order solution. The 
subscript 2 is used for the perturbation solution 
to avoid confusion with the first-order solution of 
reference 1 . 

The body shape, consistent with equations (12), 
r b -=jr m (vt>, o+e.r v ) 0<1 ( 13 a) 

r>=i n ( e ^ v ) ly(<r+I) 0=1 ( 13 b) 

where 



Equations ( 13 ) follow from equation ( 4 a) and 
the asymptotic form of equation (12b) near 
7 ] hii} . (See appendix B for further discussion.) 
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For problems wherein the body is presented, 
equations ( 13 ) are used to determine e and N. 
For these problems the constant a> in equation 
(12a) is initially unknown and is found as a 
consequence of the solution. 

Substituting equations (12) into equations ( 3 ) 
and collecting terms of order e° and e 1 yield the 
zero-order and perturbation equations, respec- 
tively, winch are summarized as follows: 

Zero-order equations: 

Continuity: 

(<Pq V ) + ^o<A) + a ^ 0 ^°= 0 (14a) 

V 


4j2 $2 | / >0 Fl\ <P> | | ( V+Um" ^2 

F 0 W() \Wo F 0 / T Uu 7? — J^n 


where 



+ + 1 ) + ft) F 2=0 
V — v^o JF 0 


( 17 c) 


1 A r 

<T + 1 


The boundary conditions at rj= 1 are (eqs. (Eo) 
and (EG)) 


— ii [1 + (r+ 1 ) m] — ( 1 ) ( 18 a) 

a 2 7+1 


^-Moment urn : 

(m-v)vo+j-(° 2 1 ')P‘*>= 0 (14b) 

Energy : 


~~~ [1 + (o"+ 1 )m] — F o(f) ( 18 b) 

a 2 7+1 

( 18 c) 

a 2 


(n~v)(p-y'Q-(°+VP=0 ( 14 c) 

where primes indicate differentiation with respect 
to 77. The boundary conditions at 77— 1 arc 

^(i)=F 0 (i)= 2 /( 7 +in 

ft(l)=( 7 +l)/(r-i) / (lo} 

Equations ( 14 ) and ( 15 ) completely define the 
zero-order flow field. The body location is found 
from the tangcncy condition (eq. ( 4 a)) which, 
for the zero-order solution, becomes 

^d(*?m) = 7 7 m ( 16 ) 

Equation (10) is, in fact, the basis for determining 
77& >0 . The zero-order solution is discussed in 
references 1 and 2. 

Perturbation equations: 

<P2 V>2 _j_ 

y—<po wo T7/77— <#)' 

^=0 Cl 7 a) 

Vo 


' ^0 + ^0+ (°’+ 1 )m 


v—<p» 


*P2 

y—<Po 


l F 0 F 2 | 
(v — <Po) 2 V+ F 0 ' 



) — <Po 

<P2 

L V~<Po J 

V <po 


+ 




( 17 b) 


where, from reference 1, 

' K = ^ y ] r ,) 2 [ 3 ( 7 + 1 ) ( ff + 1 ) P — 4 < tt ] 

^“( 1 ) = ( 7 - 1 )^ 7 + 1)^ < ' 27_1 ' 1 ( <7 + 1 )/3 

— 2(77 (7— 1)] 


^(1)=^+^, [ 3 ( 7 + 1 ) (<T+1) $ 2(7 (7- 1)] 

When a body shape is specified, the constant a 2 is 
initially unknown. Its value must be such that 
the tangent flow boundary condition at the body 
surface, namely (eqs. (Bll) and (ElO)), 

<P2{Vh,o)^(<r+l) (y+yv—P)/y / 3 <l ( 19 a) 

lim (i?'^) = ( 7 + 7 M— 1 )/t 0=1 ( 19 b) 

»j -»0 

is satisfied. 

Since the unknown constant a 2 appears in the 
boundary conditions at 77=1 (eqs. ( 18 )), it is not 
possible numerically to integrate equations ( 17 ), 
starting at 77— 1. However, if equations ( 17 ) are 
eacli divided by a 2} the quantities (p 2 la 2i F 2 /ef2, 
and ^ 2 ja 2 can be considered as the dependent 
variables, and a numerical integration of equations 
( 17 ), starting from 77^=1, is then possible (using 
the boundary conditions listed in eqs. ( 18 )). The 
numerical integration cannot proceed all the way 
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to 7? ft>0 since the differential equations are singular 
at Unit point. Inst (aid, the numerical integration 
is terminated at some point near for which 
the asymptotic solutions derived in appendix 
D are valid. If the values of FJa 2 , fa!a 2} 

y—<Po, fa, Fu, and 0 o (defined in ref. 1 and appendix 
D) are known at this point (from the numerical 
integration), equations (D7), (D12), or (D17) (de- 
pending on whether O<0<I, 0=0, or 0=1, 
respectively) ran be solved 2 for the constants a 2 , 
B, and D. (The constant a 2 is introduced into 
equations (D7), (1)12), or (Dl7) hy replacing 
<P 2 , fa } and F 2 by a 2 (<pja 2 ), a 2 (fala 2 ), and a 2 (F 2 /a 2 ), 
respectively.) When B and D are known, equa- 
tions (1)7), (Dl2), and (D 17) completely define 
the flow in the vicinity of 7j hw0 . When a 2 is 
known, the dependent variables of the numerical 
integrations ran he converted hack to <p 2 , fa, 
and F 2 and the solution is complete. 

EXPRESSIONS FOR SHOCK SHAPE, PRESSURE DISTRIBUTION, 
AND DRAG 

Expressions for shock shape, pressure distribu- 
tion, and drag are summarized in the present 
section in terms of dimensional quantities. The 
cases 0<1 and 0=1 are treated separately. 

Case 0<1, — The zero-order body is of the form 


r t >,Q(x) = rj t , t0 C7 m (20a) 

which defines C for a specified body. For a power 
law perturbation 


r b [jr) 
r bf o(7) 


1+— 

Vb,o 



(20b) 


which defines e for a specified body perturbation. 
The corresponding shock shape is 


( 21 ) 


2 If the loading form in 7? is the only term retained in the expression for 
F 2 and if the leading terms in A and 7? are the only terms retained in the 
expression for then 1 results 


m-L iW-'Y 

,oW v»,o \L/ 


jar) 


C*+n(7+*VM— 0 )/y 

<P3 , n(r) — v?o) / Fi\ 

0<1 

02 (y—0) Fa \ / 


(t+7M_ D/y 


~ T <pi , ti(v~v 0 ) F>~\ 

V L 02 ' (■> D Fo Oj J 
a 2 E? 

XJ ~( 7 -^)E 0 a2 


0<l 


for v near 17 6 0- 


where the zero-order shock shape is I\ 0 (x) = 
CP"=r M (5)/i? M - 

The local pressure coefficient, at any point, is 
(from eqs. (1), (5a), and (12c)) 


c»= 


P~~P« 




(22a) 

(22b) 

(22c) 


In appendix C it is shown that 

To (vt>) —F 0 (vb, 0) 

F> ( Vb) —F-iivb. 0 ) (23) 

so that the pressure at 77^0 is the same as that at 
7] b for a given 7. Using equations (22b) and (23) 
gives the pressure coefficient on the perturbed 
body : 

f/>,6 2U 0 ( 77 ^ 0 ) r . . F 2 (T}t )t a ) / UUH 704\ 

(dr &f0 /d.rj 2 i?&,o L \L/ J 


Alternative forms can also be deduced. 

The forehody drag up to station 7 can be found 
hy integrating the pressure distribution along the 
body surface. This is done in appendix E and 
the result is (from eq. (E4)) 

DCr) = _ 2F 0 ( Vh , 0 ) f/7 

2**m*q6 2 \r t , t2 (L)]'+ 1 ™W+3)-2\V/V 

e[/n(<r+3) — 2] fF 2 (r} bt 0 ) 
m(<r+3) — 2-j-A' L^o( , ?6 1 o) 

| f<r+ !)(!+/*) j 

Equation (25) is valid provided that 

m>2/(<r+3) 

.V>2-m((7+3) 

Equation (26a) corresponds to 0<Ch which is the 
case under consideration in the present section. 
For tr = 0, equation (25) assumes that the body is 
symmetric about the 7-axis, The overall fore- 


7 (25) 

(26a) 

(2Gb) 


527024 GO 2 
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body drag coefficient, referenced to the cross- 
sectional area of the base of the zero-order body, is 




im _ 


+1 


2^m WoQ ^o) / c[mfr+3)-2] 

m(<r+3)— 2 \ ' m(cr+3)— 2+ A r 

F 2 (^,q) (cr + 1 ) ( 1 +m) H \ ^7) 

F 0 {Vb r o) Vo,o JJ 


When the perturbations are due to boundary- 
layer development on the zero-order body, the 
pressure drag on the zero-order body can be 
found from equations (25) and (27) by omitting 
the term (<r+l)(l + m)/ 7 ?»,o* 

The drag can also be found from consideration 
of the energy of the transverse flow field (refs. 
1 and 2). This is done in appendix E. The 
resulting expressions for drag (eqs. (Ell) to (E17)) 
are more general than equations (25) and (27) 
since the former can be made to apply for all 
values of m and 2V. 

Case 0—1, — For 0—1, the zero-order flow cor- 
responds to flow over a flat plate (<r— 0) or cir- 
cular cylinder (o-=l) of semi thickness, or radius, 
equal to r x . If the nose drag at x— 0 is known, 
and is denoted by T? X) the zero-order shock shape 
is given by (ref. 1) 

J?o(5 
r x 

where C 1>N =D N !2'~ (r Tr ff q(r N y+' is the nose drag 
coefficient and I is tabulated in reference 1 as a 


! j >+3 


\ 2/ (<r+3) 


(28) 


r 



Perturbed constant-energy flow. (Flow due to drag im- 
pulse at 7=0 plus perturbations associated with 7 b {x)). 


r 



function of a and 7. The body shape can be 
expressed as (from eq. (13b)) 


r b {7)_n oQ r f r 3 Ay7 " +1) 

r v |_ \ /> r x / _ 


(29) 


which defines e and N for a specified zero-order 
flow and body perturbation. The body ordinate 
r 6 (x) is measured from the surface of the zero- 
order body (sketches (a) and (b)). Sketch (a) 
is the hypcrsonic-slendcr-body-theory idealiza- 
tion of tlie flow pictured in sketch (b). Hyper- 
sonic slender body theory gives a poor represen- 
tation of the flow near x=0, as previously noted. 

From equation (29) it is seen that e(x/L) iV = 
[r 6 (7)//?o(x)] ff+1 . The perturbed shock location 
can then be expressed as 



From equation (22c), the pressure coefficient on 
the body is 


Cp, b 

(dli u jdx) 2 


—2F u (d) 



(31) 


The forebody drag up to station x is, for A T )>0 
(from eqs. (E7) andj(28)) 


T)(x) D N 


", , (1 +AEM 

n(<r+l)I 


/ r,( x) V +1 ~ 

\/i 0 (x)J 


(32) 


The perturbation solution is valid provided 

fo(F)/^ 0 (:F)]» +, <s:]. 
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SOLUTION FOR V^O 

When A r =0, Hie solution of llie perturbation 
equations can b(* found analytically in terms of 
the zero-order flow. The procedure is as follows. 
For a shock shape given by 

ea 2 ) (33a) 

the perturbation equations and the boundary 
conditions at tin' shock are satisfied by 


<P-2 ' 


Ei 

a> 


2 Fv—nK 


(33b) 


(33c) 


J£— ntf. (33(i) 

Note that F 2 (T} t)i0 ) = 2a 2 F iy (7} bt0 ) since F f 0 (rj bt o ) = 0. 
For j3< 1 the boundary condition on the body is 
satisfied if (eq. (B9)) 

<P' 2 (Vb, o) = 1 — <Po(Vb. o) (34) 

Evaluating equation (33b) at ^ >0 , noting <a>(^o) 
~Vb,o, and substituting into equation (34) then 
give 

I/??*, o /? <: C 1 (35) 


For (3— 1, the shock perturbation 1R may be 
considered as due to a body perturbation A?y. 
From equation (28) 


&H_ <r+ 1 A7 iV 

r.v 


(36) 


But, from equation (33a), AR/R^=ATi/Ii^ea 2 . 
Substitution into equation (36) then gives 


€&2 


<r-r 1 Ar jV 
o' +3 r N 


i8=1 


(37) 


The present solution results from the fact that 
for A T =0 the perturbed shock and body follow 
the same power law as does the zero-order flow. 
The resulting flow can in fact be treated as a zero- 
order problem. However, the solution is useful 
for making a partial check on calculating machine 
programs (when the latter are used to obtain 
numerical solutions of the perturbation equations) 
and for providing additional data when tabulat- 
ing a 2 and F 2 (ij m ) as functions of p, 7, a, and 0. 


NUMERICAL RESULTS AND DISCUSSION 

The equations of motion have been integrated 
numerically to determine the effect of (1 ) boundary- 
layer development, (2) very small angles of attack 
(for (7 = 0), and (3) blunting the nose of very 
slender wedges and cones. The results arc dis- 
cussed herein. 

EFFECT OF BOUNDARY-LAYER DEVELOPMENT 

The boundary -layer displacement thickness <5* 
on a slender power law body at hypersonic speeds 
is derived in appendix F. The result shows (re- 
calling that the superscript bars have been omitted 
in eq. (Fl3)) 

3 


5* 7 — 1 / co 

mj¥y\K 


M. 


J< (l + 


J 1 4~ Ja\ 

J 2 ) ( * 


L F o( r) b 0 ) 2 rn (cr — J— 1 ) — 1 \Ij 


(38) 


where /u,Z=p 00 w 00 i//I a; , w=-TJji e ITji m = cons taut, 
and {!„, is free-stream viscosity. The quantities 
J 2 and J^)jJ 2 can be found from table I for 
specified 

q - T—l 2(1 —m) _ 7 1 2(a+\)0 

7 2m(<r+I) — 1 — 7 (crT 1 ) (4 — 0 ) — 2 


(fast) w 


(39a) 


stagnation enthalpy of fluid at wall 
stagnation enthalpy of fluid in free stream 


(39b) 


For 0<jd<l and <r=0,l, the quantity 0 varies 

between the limits 0 <£<2(7— 1)/7. Table T is 
based on the numerical results of references 6 to 
8. The assumptions involved in the derivation of 
equations (38) and (39a) are noted in appendix F. 
Equation (38) is not valid for <r“l, 0~] as dis- 
cussed in appendix F. 

The effective body shape is r A =Fa ?0 -|-5*. For 
5* small compared with r M , the effective body 
shape can be expressed as 

«»> 

where e is found from equation (38). Thus, for 
this case, 


( 41 ) 



8 


TECHNICAL REPORT R~45 — NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 


The perturbation equations have been integrated 
numerically for <r=0,l; 7— 1.15, 1.4, 1.67; 0 = 0, 
% U, %, 3 ; and values of p defined by equation 
(41). The case <7=1, 0=1 is excluded. The 
results are given in table IT. 

For m=y 4 (i.e., j<=0), the perturbed shock 
shape and body shape are similar, both following 
the % power law variation of the unperturbed flow. 
This fact was used in references 9 and 10 to study 
boundary-layer development at hypersonic speeds. 
For this case, (lie perturbation solution is given by 
equations (33). 

EFFECT OF VERY SMALL ANGLES OF ATTACK (FOR 

The effect- of very small angles of attack on two- 
dimensional power law bodies at hypersonic speeds 
can also be found. This is done herein. 

If a two-dimensional power law body is at angle 
of attack a , the equation of the upper surface 
becomes r b —7} bf oOx m — ax, or 

r b =x m (7) bt o+ a l ~ ,n ) (42) 

where 

e = — a/S 

Assuming e<Cl, the resulting flow (in the upper 
half plane) can be found from the perturbation 
equations with 

ju=0/2 (43) 

Numerical solutions have been obtained and these 
are tabulated in table II. For ju = 0=O, equations 
(33) apply. 

The lift per unit span can be found by noting 
that the perturbation solution is antisymmetric 
about the r-axis. The local lift coefficient is then 


A 

The moment about the leading edge 3f is 
.17 f x(Ac P ' h ) cl./ 

Integration yields 

= 7 -— 4 5F> (i? M ) a (46) 

I?q 1 + ™ 

The center of pressure x c p =^fjt is then given by 

L w+1 

The lift problem, for <r=l, is formulated in 
appendix G. It has also been treated in reference 
11 using Newtonian theory. 

EFFECT OF BLUNTING THE NOSE OF VERY SLENDER WEDGES 
AND CONES 

The effect of blunting the nose of very slender 
wedges and cones is now considered. It is as- 
sumed that the wedge or cone is sufficiently slender 
so that the major contribution to drag is due to 
t lie blunt nose. The zero-order flow is then a 
constant energy (0=1) flow. The divergence of 
the body downstream of the nose induces a small 
perturbation in this zero-order flow. 

Sketch (a) indicates the flow field considered, 
and sketch (b) indicates the corresponding physi- 
cal flow. In the present example, the body shape 
is given by 

r b =S b 7 (48) 


— ((*/;,?>) bottom (^p,&)top 

= — 4:€m 2 d 2 F i (r} bJi )x m ~ 1 
The net lift per unit span L is (for m 5 ^ 0 ) 

L — qL f Ae P ' b dx 

J o 

= — 4e mS 2 LqF 2 (tj br0 ) (44) 

The lift coefficient, referenced to Z, is (recalling 
€ = — a/5) 

A 

-k-— — 4e m8 2 Fi (jj 6 , 0 ) 

Lq 

(45) 


whore r b is measured from the surface of t lie zero- 
order body (sketch (b)) and 5 b is the semivertex 
angle of the wedge or cone. From equations (13b) 
and (48) it follows that 


<r+l 

2~ 


(49) 


for the present problem. The equations of motion 
have been integrated numerically for ^=(<r+ 1)/2; 
0=1; <r = 0,l; and -y = 1 .1 5, 1.4, 1.07. The 

results are given in table II. The shock shape, 
pressure distribution, and drag can be found from 
equations (30) to (32). The solution is valid 
provided that (/* 6 /J? 0 )* +1 <Cl . 


= 4w5F 2 (ij6,o)« 
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Define 1 T i e following quantities: 




(50a) 


(50b) 


(50c) 


where K is a function of y,a. Equations (30) and 
(31) can (lien be written 


— = 2 /' 0 ( 0 ) 


Si 


I! fx(7)”+ :t jnii, [ jL (7)'+3 J ' ' j- 

\&**i 


(5 Iii) 


(51b) 


The quantities II and c Ptb /8l are functions only of 
x for a given <r,y. 

The problem of a blunted wedge or cone has also 
been treated by Chernyi in references 12 and 13. 
An account of his method is given in appendix IT. 
The method is approximate, but becomes more 
exact as y approaches 1. It has the advantage 
of not being restricted to small values of J as is the 
perturbation analysis of the present report. 
Chernyi presents curves of shock shape and surface 
pressure distribution against x for 7=^ 1.4 and 
cr— 0,1. Chernyhs equations are solved analyti- 
cally for small x in appendix IT, and the resulting 
expressions for K, a 2} F 0 ( 0), and 7%(0) are compared 
in table ITT with the numerical results of the 
present report. Tin 1 agreement becomes poorer as 
y departs from 1, particularly for F 0 (0) and a 2 . 

Lewis "Research Center 

National Aeronautics and Space Administration 
Cleveland, Ottio, May 22, 1959 



APPENDIX A 

SYMBOLS 


a 2 

constant defining shock perturbation 
(eq. ( 9 )) 

C 

constant (eqs. (6)) 

t D 

drag coefficient- (cq. ( 27 )) 

c Dm 

drag coefficient, for impulsive drag ad- 
dition at 5— 0 (eq. ( 28 )) 

c p 

local pressure coefficient 

e.v 

impulsive drag addition at 5=0 

D{x) 

forebody drag up to station 5 

Fo(n),F*M 

pressure similarity variables (eq. 
(12c)) 

g 

a constant (eqs. (D6b), (Dllb), 
(DlGb)), or ratio of stagnation 
enthalpies (eq. (F 3 b)) 

(J w 

wall- to free-stream stagnation en- 
thalpy ratio (eq. ( 39 b)) 

h„ 

s t agn a t ion on tl ial py 

I 

function of <7,7 tabulated in ref. 1 


quan t i t i es defi n in g boun d ary -1 ayer 
thickness (eq. (F 9 )) 

L 

slroamwise length of body 

M 

free-stream Mach number 

m 

zero-order power law exponent (eq. 
(6)) 

N 

perturbation power law exponent de- 
fined by equation ( 9 ) 

Pr 

Fraud tl number 

V 

pressure 

q_ 

dynamic pressure, p^u*/ 2 

Rix) 

lateral coordinate of shock 

r 

lateral coordinate (normal to 5 -axis) 

r b (x) 

lateral coordinate of body 

T N 

semi thickness or radius (at- 5 — 0 ) of 
blunt-nosed two-dimensional or axi- 
symmetric body 

T 

temperature 

u 

velocity in 5-direction 

V 

velocity in z- or r-directions 

w 

velocity in 0-direction of cylindrical 
coordinate system 


Cartesian coordinates with 5 in stream 
direction and origin at nose of body 

a 

angle of attack 


alternative zero-order shock shape pa- 
rametcT, 2 (~ 1 )/(< r + 1) 


0 <'<!• (F6) 

7 ratio of specific heats _ 

<5 characteristic slope, H 0 {L)lL 

5 * boundary-layer displacement thick- 

ness 

e small quantity, defined by eqs. (20b) 

and ( 29 ) for specified body per- 
turbation 

77 lateral coordinate similarity variable, 

r/ 7 ?o (cq. (10)) 

77^0 r &# o//fo, tabulated in ref. 1 as function 

of 7, <r, and j 3 (denoted by 77 
therein) 

0 cylindrical coordinate (sketch (e), ap- 

pend i x G) 

do(rf) zero-order stream function similarity 

variable (cq. (Pi)) 

H alt ernative perturbation power law ex- 

ponent (eq. ( 17 )) 

7i „ free-stream viscosity 

1 * eq- (10)) 

p density 

a 0,1 for zero -order flows that arc two- 

dimensional or axisymmetric, re- 

spectively 

fpo(v ) similarity variables for v (eq. ( 12 b)) 
^0(77) similarity variables for p (eq. ( 12 d)) 
tiziy) similarity variables for w (eq. (G 5 )) 

Subscripts : 

b quantity evaluated at perturbed body 

surface 

b } 0 quantity evaluated at zero-order body 

surface 

s quantity evaluated at (or just down- 

stream of) perturbed shock 
0 zero-order solution 

2 perturbation quantity 

qo undisturbed free-stream value 

Superscripts: 

(-) barred quantities are dimensional 

() unbarred quantities are nondimen- 

sion al (eqs. (1)) 

()' primes indicate differentiation with 

respect to 77 
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APPENDIX B 

BOUNDARY CONDITIONS AT v = 1 AND v ^ Vb , 0 FOR PERTURBED FLOW 


It is convenient to satisfy the perturbed flow 
boundary conditions at the zero-order shock loca- 
tion (77= 1) and at the zero-order body location 
= The appropriate boundary conditions 
are found herein. 

Boundary conditions at 77= 1 : 



The zero-order shock location is at rj— 1, while 
the perturbed shock is at 77=! 4 - as indicated 
in sketch (c). The flow on the downstream side 
of the perturbed shock (denoted by subscript s) is 
given by (from eqs. ( 5 b) to ( 5 d) and ( 12 a)) 


c '- 4 t ” f '" 

. |+ "=( |+ S) {V ] 

(Bla) 

2 , . 
P~y +] {m t 


(Bib) 

T+l 

Ps y — 1 


(Blc) 


Let Q=zQ(r)) be the variation, with 77, of any 
flow quantity, at a fixed station £. Expanding in 
a Taylor series about 77=1 then gives 

QXv)=Q(\)-\-Q'(\)(ii-l)+ • - . (B2) 

If 77 is taken to be the perturbed shock location, 
1 + ea 2 £ iV > equation (B2) becomes 

(B 3 ) 


Consider Q to be the transverse velocity r. Equa- 
tions (1 2b) and (B3) then yield 

v s = [(^p 0 (l) + €{ A ^ 2 (l)] + €®2f'v , o(l)+ • * *} 

(B 4 ) 

Equating equations (Bla) and (B 4 ) and collect- 
ing terms of order e yield 



Similarly, identifying Q with p and p yields, 
respectively, 



(Bob) 

a 2 

Equations (Bo) and ( H (>) give the boundary con- 
ditions on <p 2 , F 2? and \ p 2 at 77— I and appear in the 
body of the report as equations ( 18 ). 

BOUNDARY CONDITIONS AT 71 = 71,0 

The boundary condition that the flow velocity 
be tangent to the body surface, ?’*— d/‘ 6 /df, is now 
considered. The eases and 0=1 are treated 
separately. 

Case 0<(1 . -The zero-order body is located at 
7 ^= 7 j b0j while the perturbed body is located at 
^ =T?6 0 4 -ef v as indicated in sketch (d). 



Q,= Q(\) + ea. 2 F <?'0)+ • • - 


11 
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Expanding equation (12b) in a Taylor series 
about 7 7 M , at constant £ yields 

v b =m r” 1 (^o+ { fau^o) 

+ ^o( 1 ?6,o)€f X + • « • ] + € f A ^2( J ?&,o)+ • • • } (B7) 

But for r ft =f"(i? 6i0 -l- «£ v ) 

(i+^y] (Bs) 

Equating equations (B7) and (B8) and collecting 
terms of e then give 

^a(i?6.o) = l+~ — (B9) 

For $<T (from ref. 1) 

^(i?6,o)=l” — ” (<r+l) (BIO) 

7 


and noting €f v = (e£ Ar ) 1/(<rf l) 7?£ permit to be 
written 

-r»r"HW0)+^(0)( f f) 1/fff+,) + . . .] 

+ ( e e v l 1 ^'+ 1) h;^(n»)]} (BI3) 

But 

^_„, r - , [(, + _i T D„ i , ) , , rt ,] (B , 4) 

Equating equations (Bl3) and (B14), noting 
^ 0 (Q)— 0, ^>o(0) = I/7 (from ref, 1), and collecting 
terms of order (e£ Ar ) 1/(<H_1) then yield to this order 

y*<P2(yt>) = (y+yv— l)/y (B15) 

From equat ion (D18) it is seen that equation (B15) 
can also be written 


Substitution of equation (BlO) into equation (B9) 
yields 

^.o) = («r+l)(7+W-|9)/7 (Bll) 

which appears as equation (19a) in the body of the 
report. 

Case /3 — 1. For 0— 1 , the zero-order body 
shape is r^o^O, an d the perturbed body shape is 
taken to be of the form 

^=(«r) ,/( ' +u J lB 

Expanding <p 0 in a Taylor series about 77 & t o= 0 


Urn (i?V 2) = (7+7M“1)/7 (BIG) 

17 >0 

which appears as equation (19b) in the body of the 
report. 

It was necessary to take r b of the form indicated 
by equations (Bl2) in order that equations (B13) 
and (B14) be consistent . Conversely, the required 
form for r b could be deduced from equations (B13) 
and (B14). That is, let r 6 =£ m (e£ iV )* where k is an 
unknown constant. Using this expression in 
equations (B13) and (B14) and noting lim 17^2(1?) = 

rj >0 

constant (from eq. (Dl8)) then show that these 
equations are consistent only when £=1/(0-+ 1) 



APPENDIX C 


PRESSURE DISTRIBUTION IN VICINITY OF BODY SURFACE 


The pressure distribution in the vicinity of the 
body surface is now considered. 

The local pressure is given as a function of F 0 
and Fo by equation (12c). Expanding in a 
Taylor series about 77* i0 si lows 
K(y{>)+e^F 2 (T} b ) — [Fo(r} b) o)+FQ(ri bt0 ) (r} b — rj bi0 ) 

+ • • •]+«£* [^07m) + . • •] (Cl) 


Noting / f o(T?6,o) = 0, from reference 1, and collect- 
ing terms of order e° and e l yield 


Fo(v»)=Fo(v»,o) 

F 2 (vi>)=F 2 (7} di0 ) 


(C2) 


Thus, the pressure at rj b is the same as that at 
?7*,o (for each value of £). 
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APPENDIX D 

SPECIAL INTEGRALS AND ASYMPTOTIC SOLUTIONS 


In the present section the perturbation equa- 
tions of motion are written in a form wherein the 
dependent variables are ^>2/(77— <p 0 ), F 2 /E 0 , and 
^2/Vv Some special integrals are then obtained. 
Finally, asymptotic solutions, valid near 77 — 77^0? 
are found. 

Following the procedure used to derive equations 
(55) in reference 1 allows equations (17) to be 
writ t on , respectively, 


(—)' -(t) +H—) + it f (z0-° 

V? — Vw Wo/ Oo \1 )—<Po; 00 '.Vn/ 


(Dla) 




P 1+? 

2 H- 


^_,2*oCl ( <Pi \ 

O' 

< D ‘W 

■'C”™) _ (® + <T “® I (^) + " s (»-» 

(Die) 

where 8 0 is a stream function similarity variable 
defined by (ref. 1) 

*,-«*[(.+ 1) f-fc 1 

L Jn <Po—yJ 


(Tlie zero-order quantities <p 0 , Fq, and \p Q are 
expressed in terms of 6 0 in eqs. (40) of ref. 1.) 
Integration of equations (Dla) and (Die) yields, 
respectively, 


r— f- 
Vo y — <pa J V 




-Vo or 


(D2a) 


and 


F 0 77—^ J y—<PoVo 

(D2b) 


Elimination of the indefinite integral between 
equations (D2a) and (D2b) yields the following 
special integral (of the continuity and energy 
equations) : 

pJE? -+(7/^+7— 0) j 5 — (1 + m) (constant; 0S 

V <Po ra 

(D3) 

where the constant can be evaluated in terms of the 
boundary conditions at 77— 1. If (i.e., 

A 7 =2(l — m)) and the subscript 2 is replaced by the 
subscript 1, the above equations become identical 
with the corresponding equations in reference 1. 

Asymptotic solutions of equations (Dl) valid 
near 7? « 77^0 (he., #o~0) will now be found using 
the procedure of appendix D in reference 1. 
Assume <P 2 l(y~<Po ) of the form 

-^=^[(V- A[ )Z,+ (iV+P- M )M^+ . . . ] 

77 — ^0 

(D4a) 

where A 7 , P } L 2) and M 2 arc constants. The cor- 
responding values of \p 2 and F 2 are then (from 
eqs. (D2)) 

|=^-[(.Y+ 1 )Z 2 + (tf+P+l)A/X+ • ■ .]+EM 

(D41)) 

§=0f{( 7 V+7 ~P)U 

r 0 

+ MN+P+1 )-0\MA+ . . . }+G 2 6S (D4o) 

Appropriate values of A 7 and P as well as the ratio 
M 2 IL 2 can be found by substituting equations (D4) 
into the momentum equation (eq. (Dlb)) and con- 
sidering 0 o ^O. This will be done for 0</?<l, 
/3 = 0, and /3=1, respectively. Since equations 
(Dl) arc equivalent to a third-order linear equa- 
tion, three independent asymptotic solutions can 
be found. 

CASE 0<fl<I 

In the vicinity of 0 o ~O, equation (Dlb) becomes 


where E 2 and G 2 are constants. The latter equa- 
tions correspond to equations (56) in reference 1, 
14 


d^F 2 /F 0 ) . ft 77 1 , 0 * / F2 

d#o 2Foi77 6> ol \F q \^o/ 


(D5) 
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Substitution of equations (D4) into equation (D5) 
yields tlie following three independent solutions: 


(1) For N=-(y~iy P= 1, E 2 =G 2 =Q, L 2 


~ 6 0 ' y ' [(y+ 7 M — 0 ) L 2t i 


y—<Po 


+ (/?-ym)M 2 ,A+ . . .] 


r_ 


Pa 


0 U v w [-PUi+<y+fi ^4 100 




>- (DOa) 


/'o 


0 O v T / (O+ 7 2 M 2 ., 0 o + . . . : ) 


Z 2[I 27F0 (> 7 & . 0) 


Equations (D6) can be linearly added to find the 
behavior of <^/ ( 77 — <po), fa/to, and F 2 /F 0 near 0 o ~O- 
Eeplacing the constants L 2 ,i , L 2t 2 , and E 2 by A, 
—By and D , respectively, then yields 

— =/l 0 Q ( t ) [(7+7M“^) + (0— YM)/0 O ] 

*? — Wo 

— m — ( 1 — M)#0o]+T^o(/d?o) + • * * (D7a) 

f r A 0 ; ^\-p+ (y+fif<k] + B(l -2gd 0 ) 

+Z) 0 S [1 + ( 2 + m )M o ]+ • • • (D 7 b) 

^ (y 2 fOo)Amy-P) - (2y-0) gd „] 

+D 0 £[ t( 2 +m)-/ 3 ]Mo+ - ■ • (D 7 c) 

The constant .1 can be found by satisfying equa- 
tion (Bl 1). This is done as follows. Recall from 
reference 1 that, near 0 O — 0, 


(2) For N=0,P=1 } E 2 =G 2 =0,L 2 =-L 2t2 : 

— =m4i+ Cl— m) M 2 .20 o + ... 

TJ—^0 

— L 2 2 + 22V/ 2t 2^0 + * . ■ 

Wo 

_ p^ == — (y—P)L 2 , 2+(27~/5) Af 2 , 2^0+ • • ■ 

^ 0 

d/ 2 . 2 _ (7— ft— 1 ) IV&T 7 

B‘2,2 2(27 fi) Fq (TJ&.o) 

— {7 

(D6b^ 


( 3 ) For Af= M , P=l, #2=0, Z 2 =0: 


-^^=^CO+M 2j3 0o+ . . .) 

77—^0 

^t^==^o[/? 2 ,3+(2 + /i)A/ 2f 3^0+ • - J 
Wo 

=0o { 0 + [t( 2 +m) — ^ 4,300+ • 
^ o 

4 3 a . o" 

4> 2(1+ji,)[7(^+2) — / 3 ]F 3(17^ q) 

si 


J r (D6 c) 


r / 7+ 1 7— p \ 

\\ 7-1 7 ) 

^(wJV+'-^.V) } 


7 

7-/9 


*-*>=>+1)3 

^0 


l (D 8 ) 


Evaluating equation (D7a) at ij 6i0 and applying 
equation (B 11) then yield 

7+lc+irV+ l W/ xT /Y <r /n 0 \ 

^1= t - — ” — — 2~ ’f&.o (ED) 

for the present case. 


case 0=0 

The momentum equation becomes, near 6 0 ~0 } 




+§[4+l)/i-(l+2<r)]-^- 
0 O <Pa 



(D10) 


where 


0; 


4+1) + 


7+1 

7—1 



4 
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The following three) independent solutions result: 

(1) For i?= — 1, P=2 } G 2 -^E 2 =0, 

— = 1- 




^2,1 


— 0 O 1 (0 + 2A/ 2( i^o+ - - - ) 


Zkl— a- 


F 

Mm 

L.r 


«.- , (0+2 7 M ill fl8+ • • • ) 

(7 — 1 ) [<y — )m] 

47(0-+ 1)>7|% [^2“ F u (i) 6<l ] 

= f 

(2) For #=0, P=2, E 2 =G- 2 =Q: 


] 


7+1 

7 


V (D 1 1 a) 


<P'2, 2 


n — Wo 

^2, i 


“M7>2, 2+ (2— m)M 2i 20Q+ * . • 


, £-2,2 + 3il/ 2i 2^0"f" * * * 

WO 

7^2, 2 +37jl/ 2( 2^0 4" • • • 

WO 

A/2 ,2. (7 — l)[h r- M )M 2 — 0 +2cr)/i-h<rf7 — 1 )] 

- — — 7h 


12 7 (<r+l),!r 0 [^± 1 F,(„», 0 )] 


(Dllb) 


(3) For N=n, F= 2, <7 2 =L 2 =0: 

95(0+2M;^j+ • • ■ ) 

V — Wo 

“T^ = 0O [^2+ ( 3 +/i)A/ 2>3 ^o 4 - * * • ] 

Wo 

-pr^=6 l f; [0 + 7 (3 + /i) M 2 ,3^0 + • • * ] 
iV^2,3 _ ~a(y 1) 

E 2 ~ ' f i±l 

27 (<r + 1 ) (3 + /z) ^2 F 0 (rj b ) 




(Dllc) 


Replacing the con+jants L 2tU L 2t2} and E 2 by A, B , 
and Z>, respectively, and linearly adding equations 
(Dll) then yield 

[1 + (1 — v)M] + B [— ju + (2 — m)^o1 

-}~D^o (2^o) -f - * ■ • (Dl2a) 
y=A$ 0 J (2 Jdl) +B(l+3g6l) 

Wo 

■+2»8[l + (3+M)Wg]+ . . • (Dl2b) 

y=Ado'(2-yfol)+B(y+ lygel) 

+ZJ»g[7(3+M)AflJ]+ • • • (Di2c) 

which defines the behavior of the dependent 
variables near d Q ~0. From equations (Bll), 
(D8), and (D]2a), the value of A is 

74-1 Tth-I 1 1/Y 

(o-+0(1+m) j^— g' ■"•F oOimOJ vl.o (D13) 

for the present ease. 


CASE 0 = 1 

For this case, i/ 6 ,o==0 and (for 0 a ~0) 

Oq (7— 1) 

^oV~7(<r+l) 

gfo> 7<r+1 

W ~7(«r+l) 

27+ff— 1 (7-1) (»- 1) 27+^-! 

7-1 X K 0 7 (, '+ I > '"+ 1 ) 



(D14) 


The momentum equation can then be written 

1 / [ 2 7(o , + D M— 37— 7<t+41t; 1 - j < p t 

^ ’ 2 T 2 (<r+l) n-w 


37+70-— 2 
' 2 7 z (<r +1) 



(D15) 
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The resulting three independent- asymptotic solu- 
tions are 

(» For .V— (?-‘>p_^±£ t : ) i ,a= B ,-o, 

Tj i=‘ 7^2, i j A/o— y A/-*, i ; 

{ [C-y— 1) +w]/: a . 1 

V~<Po 

+ [ty-V)+y(n-P)]\r- 2 ,i 6 Z+ . . . 


to. 


y 1 


— o 0 7 l L 2 , i — (yP~t 1 BC 1 ^ 




— e, 7 (o — 7‘/ 7 n/ 2 id' K + . . . ) 

1 0 




M 2 , I [ ( 7— I)+7m]A’ 


0 7t * +1) // 7^1 V 

!)]/\,(0j U 7 / 


Ui 7 / 7 [ 7 P— ( 7 — 1 , 

+ 9 2 / , 1 f ( 7 — 1)[27(o , + 1)m 
-7 (ff+U L 


-/ 


- 3t - t ,+ 4 I+ ^ + ^] } 


(DlOa) 


(3) For A 7 = M , F=^±^i, G 2 =L 2 =0: 


[/?,+ 0*+/’+ i)-^2,.A p + • • •] 

\K0 

4^=es {0 +[7(m-!-p+i)-i]a/*X+ • • •} " 

1 0 

- O'- D(cr-l) 

A I y 3 (37 T 7(7—2) /ve T(g t 1) 

Ey 27 2 (<t+ 1 )(Ft-m)[ 7 (m+F+ 1) — 1]F„(0) 

-A 


(Dl6c) 


Replacing 7^ 2(] , Z 2 ,2> and 7? 2 by - 1, R, and 77, respec- 
tively, gives 

<p2 =Ae 0 ( 7 ^{(7 — i+7m)+[t-“1 

V <Po 

+7^M~^ > ;]/^o } +^[~M+ (I*~ f*) ( jQu] 

+DOZ(rhOZ)+ . . . (Dl7a) 


(2) For A=0,P= 2t+ * A E 2 —G 2 =0: 

7 (cr+l) 

— a<^%2+ (7 > — m) d/ 2 , 2 0 ( 7+ • * • 


“^7 ^2,2 + (7*+ 1 ) A/ 2i2 ^+ . . . 

4*= ( 7 - - 1 ) L. 2 + [7 (C+ 1 ) - 1 ] -lA 2 .A r +. • • 

1 0 

-U 2 , 2 
Pi, 2 

(7- 1)(J - 1) 

{ju(7— 1 ) [27(0-+!) /2 — .37— (77 + 4 ]}/Vq 7(t + !> 
2 T 2 ( < r+l)7 > [7(P+l)-]]F u (()) 

(7-l)(g- ) 

[(2— 7) (37+70-— 2)] K 0 7( l f+1) 

■ t ’ 2 7 2 (<r+l)F[ 7 (C+l)-l]F o ( 0 ) 

= 9 

(DlGb) 


^=Ae 0 ( r y ‘ ) [- 1 - ( T P + 1 )ffl] + 7? [1 + (7»+ 1 ) 

Y0 

+ 770jJ[l + ( M +r+1)/<H- . • . (Dl7b) 
Wo ^ ^ ( -y 2 PM) + B{ ( 7 - 1 ) + [7 (P+ 1 ) 

* 0 

-i].<}+/^s[7(m+/M + . . . cr>i7<-) 

The constant .1 is found as follows. Considering 
only the leading terms, equation (D17a) can be ex- 
pressed as 

< D1 «) 

7 

Noting equation (BIO), it is seen that 

^7 i 7 ± -![ 2 T i ^ ( °^r (Di9) 



APPENDIX E 

DERIVATION OF EXPRESSIONS FOR DRAG 


Tlu* drag on the portion of the body upstream 
of any section x , designated 77(7), can be found 
either by integrating the pressure distribution 
along the body or by computing the energy of 
the transverse flow at section 7. Both methods 
are used herein to obtain expressions for 77(7). 


INTEGRATION OF PRESSURES ALONG BODY SURFACE 

The cases 0<[ 1 and 0=1 are treated separately. 
Case 0<T. The forebodv drag is given by 


D(x) 


r 'b& 

=27r ff Cp' b (r b )*dr b 

J o 


(El) 


assuming, for <7=0, that the body is symmetrical 
about the 7-axis. But, from equations (13a), 
(22a), and (23), 

7 b —8Lx m (i? 6(0 +€X iV ) 

dT b ^dLm x m ~ l ^,o+c J V J dx 

[FofaJ + tFiiriw) <*] 

Substituting into equation (El) gives 


D(7) 


2rr a m 3 qb^ 3 (Lrjh^Y 


+ i 


= 27^o (vt>,o) 


j: 


(a -|-3) — 3 


Fiivi,. o) j 

<r+l+— 

b W 

,.w(<r+3)-3-f A r 1 

_Ao (ijil.o) 




Assuming 


m>2/(<r+3) (i.o., /3<1) 
A T >2— w(<r+3) 


(E3a) 

(E3b) 


integration of equation (E2) yields 

D(z) 

27r' r m 3 ^6" +3 (7^ij ( ,,o)" +1 


r 


_2F«(v„, o) J 
m(«r+3)-2l 


jm f>+3) 


t[m (>+3) — 2] 
m (<r+3)— 2+ A r 


X 


_Ao(ij(),o) 


<r+l 


N 

in 

Vb, 0 _J 


„77l (V-j-3) — 2-j-A 


y (FA) 


j 


Case 0=1. For 0=1, the zero-order flow field 
corresponds to the flow over a flat plate (<r=0) or 
circular cylinder (<r=l) with semithickness or 
radius, respectively, equal to r jV . Designate t lie 
known zero-order drag of the nose (x= 0) by 77 v 
as in reference 1 . Then 

D(x) =D y +2ir a q c v ^/ b ydr b (E5) 
Jo 

where 7 b is measured from the surface of the zero- 
order body. From equations (13b) and (22a), 


“ <T+i v . m+ ff+i 


r b ~8Le (T+l x 


(E6a) 


1 m AT l 

dr b = m (l dx 


c Pib = 28 2 m 2 x 2 <m ~~ 1} F 0 (0) + 0 (*) (EGe) 

Substituting into equation (Eo) and integrating 
yield, for A 7 >0, 


D{x) D N 


lG7r^(l+ M )5 2 F 0 (0)[7/ 0 (i)] ff+1 6 


(<7+3) 2 (<7-f 1 )m 
wl i or ein m — 2 / (<r + 3) . 


(E7) 
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DRAG FROM CONSIDERATION OF ENERGY OF TRANSVERSE 
FLOW 

The forebody drag T)(x) can also be found hy 
computing the energy of the transverse flow field 
as discussed in references 1 and 2. 

Let E be the energy (computed as a perturbation 
from t lie undisturbed frec-stream value) per unit 
mass at any point. Then 


pE=p[c r (T-T J+|] 

<r’.G-fcK] 


From references 1 and 2 


T)(f) ~2ir a I pRv* 


(m 


wliere the integration is conducted at station x. 
But 


J=SijX m 7j 

d7 =dT,i m dr? 

Vs—l+ea 2 r v 

=( eJ *y 


P<i 


Substituting equations (12) and (FlO) into equa- 
tion (E9) and neglecting higlier order terms yield 


T) (7:) 


4 <pr 0 b 2 m 2 (6 L) ° + V (<J f 3) m 2 


-jrww 


-re 

^ VO v 


f F 0 , ^0 <Po 

^b, 0 


•,-i+T)"'' 1 ’ (eii) 


But, from equations (24) and (29) of reference 1 


V. + +o) r<H 



/3<1 


,3=1 


y (Ki2) 


where I is tabulated as a function of y and a in 
reference 1. From the mean value theorem and 
equations (15), 


Similarly, 


^o(O) 


(<7+1) (7-1) 


=«a 2 J' v l 


\7 1 




~ 7 S — 1 

- e.r' v 

0<+ 

1) ‘ ,,A 

0-1 

J 


(El 3) 


r (EH) 


Using the asymptotic forms noted in equations 
(43) of reference 1 it can be shown 

r, f b imv" <l»?=4 +| vVu 2 

“»/ Vb , 0 * ' 1 

|+ 1 /o(,,o)] « r ' v for ft— 0 ” 


=lii<rh('i - order 


Define 


for /3>0 J 
(El 5) 


-L [A+^(l+«)] ( E,c » 

Substituting equations (F12) to (ElG) into equa- 
tion (Ell) permits the latter to be written for 
£ = 0 ; 0<+/3<O, an <l 0=1, respectively. 


42rr <T 5 2 m 2 (5Z) <7+1 i‘ (ff+3)m '“ 2 




.V^o(^o) 


+ ei 


Vv 


Vb'oR [)(yb, o) 


<7+1 ■ ’ 1 y £ — l 7 - 

+ [ 7 + EoU.o)]’ /7 1 for /3— 0 

(El 7a) 
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_ 7 ^J}bj^o(Vb,o) I -v J” j _i_ 4a 2 

” m" (a +3)— 2 6,1 L 7-1 



for 0 <j 3<1 (El7b) 

FM 1 

(a+l)(7-l)J 

for /?=1 (El 7c) 


All the quantities on the right side of equations 
(17) are known from the zero -order solution except 
a 2 and which are found from the perturbation 
solution. 


From equations (El 7) it is seen that, when 
m = 2/(cr + 3) (i.e., /3=1), the zero-order drag rises 
discontinuously at ar~0 and is constant for x^Q. 
This flow was referred to as a “constant energy” 
flow in reference 1. Similarly, when N—2 — m 
(<r+3), the portion of tin 4 drag associated with the 
perturbation solution rises discontinuously at y— 0 
and is constant for x>0. Equations (E17) are 
more general than equations (E4) and (E7), since 
the range of N is not restricted in the former 
equations. 


APPENDIX F 


BOUNDARY-LAYER GROWTH ON SLENDER POWER LAW BODIES AT HYPERSONIC SPEEDS 


The boundary-layer displacement thickness on 
slender power law bodies at hypersonic speeds 
was found in reference 2 using the local similarity 
concept of reference 5. It was assumed in refer- 
ence 2 that <r =- 1, Prandtl number = 1, the body 
is insulated, and viscosity is proportional to tem- 
perature. A similar derivation for displacement 
thickness is presented herein excepl that the a — 0 
case is included and the assumption of an insu- 
lated body is replaced by the less restrictive 
assumption that the body surface temperature is 
the same at all stations along the body. All 
physical quantities are dimensional in the present 
section, and the superscript bar is omitted for 
convenience. In the present section the symbol 
p. represents viscosity. 

Following reference 5, the independent variables 
5 and y (where s is distance along the body and y 
is distance normal to the body) are replaced by 


PfU r 

v (2 *J 1/2 



Ay 


(Fla) 


where h st is the local stagnation enthalpy. Tf it 
is further assumed that 

Fr= 1 (F4a) 


PePc 


(F4b) 


the boundary-layer momentum and energy equa- 
tions become, respectively, 


and 


<r+j<t'- o 


with boundary conditions 


(Fo) 


/(0)-/'(0)^0;/ , (^) >1 as rj >co 

and 

< 7(0 )=g u . or < 7 '( 0 ) = 0j g(rj) >1 as r\ ^oo 


Here 




2 s (Uf e 

XL d s 


(FO) 



S = \ JlcPeVA* 0 <l» 
t/o 


(F 1 b) 

The subscript c 

represents local frcc-stream con- 

ditions just outside the boundary layer. 

A stream 

function exists such that 




C)\p 

= 0 
C )'/ 

► 

(F2) 






pl> "' 0 

1 


Assume 





t=(2syy(i) 


(F3a) 

so that f' = uju e 

. Similarly, assume 



7 

("S t) € 


(F:jl>) 


Equations (Fo) were solved numerically in refer- 
ences G and 7 for various values of y( 0) and (3. 
[Note: S=g—\ therein.] The displacement thick- 
ness for these solutions is (ref. 5) 



but 

l -I ^ 1 a/*) a/jcd* (F8) 


Substituting equation (FS) into equation (F?) 
permits the latter to he written 


=-^ 2 ro+F 1 u 0 (j,+ ? ahj ^ 

rlj, p, y , Lv z / - J 


(FO) 
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where 


J, 


r (1-f'w 

Jo 

2- rro-rw 

Jo 

,-J. Of— i)<ln 


J. 


For a given </(()) and Ji the integrals ._/> an<l (-/, f- 
r/. ( ) can be evaluated from the numerical results 
in tables IT of references 6 and 7, In particular, 


permit equation (Fib) to be integrated. Sub- 
stitution into equation (F 7 ) then yields, for 
(t — 1 )M *»1 and (y-l)il/*» 1 , 


a* 7-j [JT ^_( 1 

L ^ 27 A Kc,l 8^F 0 (r) bfi) ) m\2m (cr+ 1) — 1 \L/ 

(FI 3) 

whore Similarly, since 

V e 1 


Me- 


a* </ e a. 



These results are summarized in table T of the 
presen t report , Add it ional values for the insulated 
wall case — are also included in table T, 
as obtained from table I of reference 8 . 

The expression for 5* and 0 can be simplified 
foj* the case of hypersonic flow over slender power 
law bodies. For such bodies 

r b. 0 = Vb, oM j 
Pe -1^2 n) / d / 7 , n \ 2 

Too % ” vlo \ d xj J 

Consistent 1 with equation (F 4 b), a constant oj 
can be defined such that tx/T=ufx m jT m or 


Y (FID 



Moo 

T„ 


(F12) 


(If fi c is not proporl ional to T t) a mean value of 
w is used.) Substituting equations (Fll) and 
(F 12 ) into equation (Fib) and assuming $=x 


§^ f 2r» ( ff +V-l 


dff e 7 — 1 d (]hlpj) (isentropic flow external 
a e 2y (pelpa.) to boundary layer) 

equation (FO) becomes 


q___ t 1 2(1 — m) 

7 2w(cr-fl)— 1 


(F 1 4) 


Equations (F13) and (Fl4), together with table 
T, define the boundary -layer development on a 
slender power law body at hypersonic speeds. 
There appear to be several misprints in the final 
expression for 5*/Z in reference 2 . (The integrals 
J v and J 2 are missing and m is in the numerator 
instead of the denominator therein.) 

The validity of equations (FI 3 ) may be ques- 
tioned for 0—\, <r^=l, since r M — 0 for this case 
and the rig] it sides of equations (Fl) and (F 2 ) 
are identically zero. If r ft<0 is replaced by r x in 
equations (Fl) and (F 2 ) and if the assumptions 
incorporated in equations (Fll) and (F 12 ) arc 
again applied to permit integration of equation 
(FIG), it is found that * lias a logarithmic singu- 
larity at Z“ 0 . It thus appears that the hyper- 
sonic slender body approximations must be aban- 
doned in order to obtain a boundary -layer 
solution for the 0~\ } a=\ ease. Such a develop- 
ment is beyond the scope of the present report. 


APPENDIX G 

AXIALLY SYMMETRIC BODIES AT VERY SMALL ANGLES OF ATTACK 


In (lie present section equations are developed 
which permit finding the effect of a power law 
lateral perturbation of t lie centerline of a slender 
power law body in a hypersonic stream. The 
flow about axially symmetric power law bodies 
at very small angles of attack is probably the 
most important application of these equations. 
The case a— 1 is excluded from consideration. 

First, the equations of motion will be put in 
cylindrical coordinate form. Let ( x } r } 6 ) be a 


Energy: 



d 

dr 



-0 


(Gld) 


Equations (Gl) reduce to equations (3) for the 
case of two-dimensional or axisymmctric flow. 

Assume a zero-order body and shock of the 
form r btQ — 7 ) b>Q x m and F 0 =x m , respectively. If the 
centerline of the body is displaced locally by a 
small amount. 

2=e.r A ' +m (G2) 



in the x,z plane, the new body and shock locations 
are given by 

r b ^x m {rj tf ^+€jr y sin 0) (G3a) 

J?= x ni ( 1 + ea,r x sin 6) (G3b) 


where a 2 is a constant. For the case of a body at 
positive angles of attack a, 

— a/8 (G4a) 


N=\-m m \ 

p—fi/2 J 


(G4b) 


cylindrical coordinate system such that x is in 
the free-stream direction (see sketch (e)). Let 
v,w be the local velocities in the r,0 directions, 
respectively. Dependent and independent vari- 
ables are assumed to be nondimcnsional (e.g., eqs. 
(1) and (7) of the present report and eqs. (8a) 
(8b) of ref. 4). The hypersonic slender body 
equations of motion (e.g., eqs. (9) and (1.0) of 
ref. 4) become, in cylindrical coordinates, 
Continuity: 

dp dpr opr + !^ =0 (Gla) 

dx~ dr r T r 50 


/-Momentum : 

dr. ^ i ^ 

dx 1 dr' r dO r p dr 


(Gib) 


0-Moment u in : 

dw dw w dw vw . 1 d//_ 
dx' 1 dr~*~ r dd r pr dO 


(Glc) 


From equation (G4a) it is seen that a must bo 
small compared with 8. 

Equations (Gl) can be reduced to ordinary 
differential equations as follows. Introduce new 
independent variables £,tj as defined by equations 
(10). The dependent variables may be expressed 


as 

v= m£ m l (<A)+e£ A shi 6 <p 2 ) "1 


w=mef v+w - 1 cos 0 0 2 
/j=(mf , "“ 1 ) 2 (/ r o+ef v sin 6 Fo) 
P~'l / 0 


(Go) 


where <p f F, yp, and Q are functions of 77. Equations 
(G5) are consistent with the boundary conditions 
at the shock and at the body surface. Substitution 
of equations (G5) into equations (Gl) gives zero- 
order equations identical with equations (14) and 
perturbation equations of the form 
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Continuity 
<P 2 ^2 


Wn n / 71 — <£ 


?7 — tAo Wo V ■ V — <Po 


+ 


<Po~\ <Po~{- (<7 + 


^ — <A) 

/-Momentum: 


^0 


0, 


^|fi| — <Po) 


=0 (G6a) 


<^2 


1 


_FoF± 

V <Pv (v Vo)‘ +) ^0 

(ff+1) 


+ 


0-Momentum: 


(2 -•*)- 

<Po 

r f \ ife 

L^-wVoJ fa 


V <Po 

+ 


<Pi 

V-fo 


N =0 (GGh) 


«2 + 


(<r+l) 


»(H- 


77 < p Q 


F, 


vfa(v~<Po) Fo 


, — 0 


Energy: 


'2 1^2 1 / *o\ *>2 . p 

p 7 7~+ ( 7 7 77 ) ~ r7 

0 ro \ yu 7- q/ 77 — ^>0 jj 




| ( o~+1)m ~ fa 
'/'(I V—<Po J &> 


$> 1 (<r 4 - 1 ) 0 *— j 3 )*l F-i n , r .. n 
— 77-+ 77=0 (GOd) 

_ ro V v^o J ' 0 


The boundary conditions at the shock can be 
shown to be satisfied if 


1*0) 


<*2 
f 2 (1) 


7+1 


fl j (a + l)/d“WoO) 


fl 3 — y+\ t l + ( <r + 1 )M]-^( 1 ) 


fc(i) 




a,(D -2 

7+1 


a* 


(G 7 n) 

(G 7 b) 

(G7c) 

(G 7 «l) 


The tangent flow boundary condition at the body 
surface is satisfied if 

*P‘2 0 ? b, 0) — (<r + 1 ) (t+W— P)Jy / 3 < 1 (G8a) 

^?o( 0 ) = (t+tm — 0 /y o'—O, 0=1 (G8b) 

It can also be shown that F 2 (tji,) — /^(tjm)- 


Asymptotic expressions for the dependent var- 
iables, valid near 170,0, are required if equations 
(GC) arc to be numerically integrated. The der- 
ivation of these asymptotic expressions is now 
outlined. The cases 0 <j 3<1 and 0^0 are con- 
sidered. (The case o=0, > 3=1 can be found from 
the equations in the body of the report, and the 
case o=l, j 3 =l is not being considered.) First, 
the equations of motion are written in the form 
(treating tpzii 7 } — <pa)t FJF 0, v+'+u, and fL + as the 
dependent variables and 0u as the independent 
variable) : 

Continuity: 

dlW(iy— gp] _ 

d0 o d0 o 

+T+ _+M +riW 0 (G9a) 

OoL*l — <Pv Yv > (ff+U^J 


•/-Momentum : 

d[<A?/(iy— ^o)l ■ k\ <p> h d(F s jF 0 ) 

d0o ^00 0o " f/3/7> d0„ 


_i_ 0 

+ 05 - , ' ,1 "Vrt> /v 


for 0</3<1 


*Po)] , N <p-> h A(FJF 0 ) 


<100 


'+ 


do V -<Po d\ lift, 


+ 


fi' l ( fa 7 * ;\ 

8 \fa Fj 


©-Momentum: 

(l(O^) , A',O a 


0 for (8 -0 




(G 9 b) 


<100 00 V 0j +WT) L 


(G9c) 


Energy: 

d(fyn) <i(^o) . 1 ( 

Ado 7 d0b + ft 


1 (a Vi ^2 Fj\ 

.1/3" -+7/<7 — M 1= 

'o\ V — <Po fa I 0/ 


0 

(G9d) 


where, for 0 o ~O and 0 </ 3 <( 1 , the quantities 
. . A’ 5 are constants given by 

*.=§-+^§-+<7)-^ ««« 


Z±1 

T 


'AN 0 ’ S " S "'™ -"f [ N 1 '■'«(’•■*>] 

(GlOb) 
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y (Gioc) 


i K eJlz . (7-1) foCnt.o) r 2 

ff+1 v 2 to (<i+l)(7 + l) 1 ?§.iiL(7+l)^ r o(’j6,ii). 

(GlOd) 


~ |i/> 

)J 


l-2(^) 1_ 

As 2 <r+l 2 M <r-r1 


(GlOe) 


case o:/?<! 

Assume (lie dopendonl variables lo beof the form 


/<; 


<P2 


y—<Po 


?=^+ l [F m +^-vmF n +...) 


=6'o [^jo A 0l sn ' l ‘Pn~\~ • • • ] 


y (Gn) 


-*=«na*o+C v/1r> n„+ . . .] 

V 

f=eh*2o+e l 0 -v m t 2 i+ • • •] 

where / yT 2o, F 2i , <P 2 o, <P 2 i, and so forth, arc constants. 
Substituting equations (Gil) into equations (G9) 
and collecting the coefficients of the lowest order 
tonus, and the second lowest order terms, yield, 
respectively, 


(A -f- 1)^20 (A m)^20 


S2 2 


<r+l 

~ k-2 ( A ~r 1 ) F 20 d - ^*3^20 — 0 

(A 7 -\-k$) O20 ~0 
7 (fJ. — A 7 )&>0 + P<p20= 0 


=0 


(Gl2a) 


(P 1 2 ~y) 1 -y-^) ^‘-^+r =0 

(5' :-k'i)tp>o — C f A +2 — F Jl -i~k^f'n=0 

(^+l-f+t*) ^.-^0=0 

7 (/i+f- A- 1 ) fci + tel =0 


(G 1 2b) 


From equations (Gl2a), a nontrivial solution 
results if 

(Cl 3) 


These define four independent asymptotic solu- 
tions. The constant coefficients in equations 
(Gl 1) referenced to F 20 may be expressed as 


^20 frgOV-r i ) 

F>20 k 3 

<P2o _—y(ii—N) \f/ 2 o 


Fn 


P F 2 < 


jr= (*+ 1 ) r ( N+ 1 ) jF- (V- M) >-°l 

I 1 20 L ‘ 20 1 20 J 


ft. 


F i0 ~r , , P 


N+\-^+h 

7 


W21 


S2 ai 


/,2 ° (a+ i)(y+2-^j ^ 20 




WlM 

20 


F. 2] 


(*•• > 0 F 


(G14) 


/ '’ 20 /%^A r +2-^ 


CASE 0 = 0 

Assume the dependent variables to be of the form 
Fo 


F 0 


=er 2 (F 2 0 +elF. u +. . .) 


<Pl 


V Wo 


-^((^20 + ^0^21+- - •) 


a 


h (G15) 


— #0 (^20+ 00^21 "T • * •) 
V 

® o (^20 + #0^2 1 + . • . ) 
Wo 


Substituting into equations (G9) and collecting 
the coefficients of the lowest order terms, and the 
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second lowest order terms, yield, respectively, 


(A T h 1 )*»- =0 

(J-f- 1 

( A I k\)<Px>~ ( A 7 T2)F 20 + i*3^20 — 0 
(A +£ 5 ) 

Cm A ) ^20 ~ 0 ^ 


(GlGa) 


(A"+ 3)^2i— (-V+2— ju)^ 21 — 4” ^=0 

cr^ 1 

(A + 2 + £j)y?21 — A 4) F 21 + ^3^21 — ^3^20 — ^ ». 

(A 4“ 2 -f" k$) f22i AqF>o= 0 

(AM- 2 — l i)( F 2 q — yh\) = 0 , 

(G 1 6 b) 


From equations (GlOa), a nontrivial solution 
exists if 


(G 17 ) 


These define four independent asymptotic solu- 
tions. The coefficients in equations (G 1 5 ) , refer- 


enced to F 20 , are 
<P 2 o h(N+ 2 ) 




A +A7 


^0 


for N=-h, - 
for iV - fx , —2 


hn A 1 2 (A -f 2 ) 

F 20 

-0 


for A 7 =a* 
for N= — k- 


1,-2 


1)(A 7 -; l)(A 7 +2)/,u 

F20” 


-0 


AM /■, 
F, 


for X=-k, 


for N--n, — 1 , — 2 


A’ 4 


20 .V+2 + A-: 


^21 


20 a 


1 r 

^ — (; 

V +3 L 


Jfcl = l 
F 2 o 7 


A-i 2-/1) :: 


hi _ 

"F,n 


1 12 , 


0*+ 1 F, 


i,r 

r/ 

20 _ 


r (V i 2-rl*,) fi+t, jP- — A’a”] 
7 -’0 /-.(A : 4 ) L /'20 /'jo JJ 


(G 18 


APPENDIX H 

APPROXIMATE SOLUTION FOR HYPERSONIC FLOW OVER SLENDER BLUNTED WEDGES AND CONES 


Chernyi, in references 12 and 13, has developed 
an approximate solution for hypersonic flow over 
slender blunted wedges and cones. For the case 
of very slender wedges and cones, his approximate 
equations can be integrated in closed form and 
the results can be compared with those from the 
zero-order and perturbation equations of the 
present report. This is done herein. 

First the derivation of Chernyhs equations will 
be outlined. The forebody drag up to station x 
is equal to the energy of the transverse flow at that 
station. Thus 

(HI) 

where D N is the impulsive drag at 7—0. For the 
present problem 

7*=5 6 7 (FT2) 

Introduce the following nondimensional quantities: 


the disturbed flow is concentrated in the shock 
wave and moves with it. This assumption is 
nearly valid for strong shocks provided y is near 1 . 
Thus, Chernyi takes 


7=d/?/d7 
= /? 


(Thin) 


and (from conservation of momentum of the mass 
concentrated in the shock wave) 


^=— L- -1 ux-'ii) 
(2 ny d * 


(IlGb) 


Substituting equations (I1G) into equations (TT5) 
yields 




jp +i -x£_ d {l ,^ ir)+ iyyr (IT7) 


=_/L^ Y Tl = /2$V+ l "** 

U7+ “ 


x= 


=_ / Ml V +l r 
1 V’oJ r K 


P 


V 


yMl5lp a 

where 


p=^- 

P« 


=__ v 


(TT3) 


(7-1) (270' 

Chernyi numerically integrated equation (II 7 ) for 
7 = 1 . 4 , (T— ~ 0, I ; the results are plotted in references 
12 and 13. 

Equation (IT 7) can he solved by series expansion, 
for small x , and the results compared with those 
from the more exact zero-order and perturbation 

solutions of the present report. Thus assume t! 
of the form (for small x) 


tions (Hi) yields 



— -LT 

/ 1 

“1 


and (H3) into eqna- 

AWy7-H 


i +-j 

(ITS) 


where K and a 2 

are unknown 

constants. 

Sub- 


i r z == = i c R = i r /? :— — = 

H„ ”^ ,rx ~y 1 Jr. ’" MT+ i k 


(H4) 

Assuming p and r are constant in the integrals on 
the right side of equation (H4) and taking p=Pfy 

yield ^noting 7 6 --7, and, from continuity, 


stituting equation (A8) into equation (H7) and 
collecting terms of order x° and j++i )2/ (<f+ 3) gi vc 


K- 


err 

:?t-2 


pr° 


<l 2 


dr= H’+'/ia-l 1 


>] 


i , r== ,= Pi, ( R°* x -r+'\ , v* pi '° +l \ ,tt-n 

J 0 \ ~^—) + 2 ( h i ) (Uo) 

Chorny i further assumes that the entire mass of 


2 ( 27 + 1 ) 

27—1 


> a—0 


(HOa) 


«2 = 


8(7+1) 


< 7 — 1 
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The pressure coefficient on the body is given by 

(HlO) 

K R° dx\ / 


Substituting equations (H8) into equation (HlO) 
and comparing with equation (51b) show 


1 * o(0) — 1/2 

CT = 0,1^ 


(j — 0 

► 

12y— 1 

“8 7+1 

<7= 1 

J 


(nn) 


These expressions for I\ : a 2 , Fo(0), and F 2 {0) may 
be compared with equation (50c) and the numeri- 
cal results of table IT. Such a comparison is 
given in table ITT. Ohernyi’s approximate solution 
becomes less accurate' as y departs from 1. 

Although the method of references 12 and 13 
is inexact, particularly for y not near 1, it has the 
advantage of not being restricted to small values 
ol “ (as is the perturbation solution of the present 
repo rt). 
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TABLE T. BOUNDARY-LAYER PARAMETERS 


0 W 

P 

■h 

•/ 1 4" 
,/ 2 

0 

0. 14 

0. 405 

0. 134 


0 

. 468 

0 


. 50 

. 4235 

257 


2. 00 

. 383 

— . 538 

0. 2 

— 0. 14 

0. 504 

0. 602 


0 

. 468 

. 518 


. 50 

. 400 

. 100 


1. 50 

. 366 

. 083 

0. 6 

0. 20 

0. 554 

2. 034 


0 

. 408 

1. 555 


. 50 

. 380 

1. 185 


2. 00 

. 204 

. 750 

1 

0 

0. 470 

2. 50 


. 05 j 

. 452 

2. 53 


. 10 

. 435 

1 2. 48 

# 

. 20 

. 408 

2. 41 


. 40 

. 370 

2. 33 


. 60 

. 336 

2. 27 


. 80 

. 312 

2. 24 

2 

-0. 10 

0. 543 

5. 67 


0 

. 468 

5. 18 


. 30 

. 333 

5. 40 


. 50 

. 274 

6. 01 


1. 00 

. 1765 

7. 85 


Results for g w are from refs. 6 and 7 wherein g u .= 

P=P, J,= l>fl } H (nc . Results for g w 1 arc 

from ref. 8 wherein « J ee q } (J r \-J. <i )IJ 2 = IT, 
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TABLE II. — RESULTS OF NUMERICAL INTEGRATION OF PERTURBATION EQUATIONS 


7 



a = 

= 0 



<T - 

-l 


Case i 


a 2 

Fiivb.v 

Case 

iL 


Fzivb, h) 

1. 15 


(■) 

_ 14 

1. 10 

0. 798 

(■) 

-% 

1. 05 

0. 803 



CO 





HI 

1. 04 

1. 90 


H 

w 

M 

1. 29 

1. 18 

(») 

0 

1. 00 

I. 04 




0 

1. 12 

1. 71 







CO 

M 

1. 04 







}i 

(») 

-% 

1. 24 

1. 39 



1. 08 

1. 49 




0 

1. 17 

1. 58 

(■) 

Ms 

1. 02 

1. 99 



CO 

M 

. 974 

2. 03 






?4 


1 

1. 40 

1. 53 


0 

1. IS 

1. 31 



O) 

Mn 

1. Jo 

I. 45 

(•) 

Hu 

. 734 

1. 90 



CO 


. 093 

1. 5G 






1 

(•) 

M 

0. 325 

0. 500 

CO 

1 

0. 0518 

HB 



C' *) 


. 188 

. 483 





1. 4 


(■) 

-H 

1. 44 

0. 8G9 

(*) 

-M 

1. 13 

0. 843 



CO 


1. 20 



0 

I. 09 

1. 91 


M 

(■) 

-K 

1. 77 

1. 39 

(■) 

0 

1. 14 

1. 59 




0 

1. 32 

J. 70 







CO 

!e 

1. 14 







% 

w 

-H 

1 . 80 

1. 04 


0 

1. 19 

1. 45 





1. 44 

1. 07 

(■) 

% 


1. 87 



CO 

H 

1 . 02 

1. 90 






% 


0 

1 . 95 

1. 77 


m 

1. 38 

1. 29 



(») 

Ho 

1. 44 

1. 53 

(■) 

H n 

. 055 

1. 05 



( b ) 

H 

. 000 

1. 30 






1 

(») 

H 

0. 443 

0 . 01 1 

CO 

1 

0. 0389 




( b -) 

H 

. 217 

. 555 





1 . 07 


(") 

_ 1 1 

1. 74 

0. 934 

(*) 

-M 

l. 20 




CO 

0 

1. 33 




l. 15 

1. 80 


H 

(■) 

-H 

2. 28 

1. 55 

(■) 


1 . 22 

1. 55 





1. 52 

1. 78 







( b ) 

% 

1. 25 

1. 99 






H 

(0 

-'A 

2. 40 

1 . 82 


0 

1. 29 

1. 40 




0 

l. 71 

1. 73 

(•) 

M 

1 . 1 1 

1. 80 



CO 

H 

1 . 10 

I. 84 






3 4 



2. 48 

1. 91 


0 

1 . 53 

1. 24 




Mo 

1. 74 

1. 59 

(■) 


, 045 

1. 53 



( b ) 

% 

. 701 

1. 31 






1 

(•) 

Va 

0. 540 

0. 070 

(•) 

i 

0. 0294 

0. 828 



(h <•) 

% 

. 244 

IHHI 






a Effect of boundary -layer development. 
h Effect of .small angles of attack. 

P Effect of wedge and cone nose blunting. 
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TABLE ITT.— HYPERSONIC FLOW OVER VERY SLENDER BLUNT-NOSED WEDGES AND CONES 

(0=1, m=(*+1)/2) 




Present numerical solution 

Ohernvi approximation (eqs. (TT9) 

and (TUI)) 

<r 

7 










K 

a* 

Fo(0) 

FA 0) 

K 

«2 

FAQ) 

FAQ) > 

i 

0 

1. 15 

0. 938 

0. 188 

0. 415 

0. 483 

0. 837 

0. 220 

0. 500 

0. 549 

1. 40 

1. 23 

. 217 

. 325 

. 555 

1. 09 

. 289 

. 500 

. 724 


1. 67 

1. 49 

. 244 

. 264 

. 605 

1. 22 

. 347 

. 500 

. 866 

1 

1. 15 

l. 04 

0. 0518 

0. 411 

0. 503 

1. 01 

0. 0756 

0. 500 

0. 529 


1. 40 

1. 34 

. 0389 

! .311 

. 643 

1. 23 

, 0938 

. 500 ] 

j .656 [ 


1. 67 

1. 54 

. 0294 | 

. 241 

. 828 

1. 34 

. 110 

. 500 1 

1 

.764 
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